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A stochastic one-dimensional cellular automaton with long range spatial interactions is intro- 
duced. In this model the probability state at time t of a given site depends on the state of all the 
^ ' other sites at time t — 1 through a power law of the type 1/r", r being the distance between sites. 

O . For a ^ oo this model reduces to the Domany-Kinzel cellular automaton. The dynamical phase 

■ diagram is analyzed using Monte Carlo simulations for < a < oo. We found the existence of two 
different regimes: one for < a < 1 and the other for a > 1. It is shown that in the first regime 
the phase diagram becomes independent of a. Regarding the frozen-active phase transition in this 
regime, a strong evidence is found that the mean field prediction for this model becomes exact, a 
result already encountered in magnetic systems. It is also shown that, for replicas evolving under 

■ the same noise, the long-range interactions fully suppress the spreading of damage for < a < 1. 
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I. INTRODUCTION 



■ The Doniany-Kinzel cellular automatorE (DK herein) has been vastly studied in the last years. Its combined 
simplicity and rich dynamical behavior make it a prototype for studying and testing different theoretical aspects of 

o 



I complex systems, such as dynamical critical propertiesH, directed percolatioroj and spreading of damageoB properties 
The DK model consists in a one dimensional lattice with binary variables associated to each site i {Si = 0: "dry" or 
"death" state and Si — 1: "active" state). The system evolves at discrete time, with the state of each site depending 
^ • stochastically on the configuration of its two nearest neighbor sites at the previous time, independent of the other 
sites at the same time. In other words, the interactions in the DK model are local or short ranged both in space and 
time. Then, the dynamics of the model is defined by four independent parameters, namely, the transition probabilities 
corresponding to the four different configurations of the nearest neighbors of any site. 
. A simple version of the model is obtained by equaling to zero the probability of any site to be in the active 
' state given that its neighbors were both in the death state (which defines a so-called "legal" rule of evolution), and 
CNJ , considering only one value for the probability of the site to be in the active state given that only one of its neighbors 
■ were in the active state, independent of which one. This leads to a two-parameter model which contains most of the 
complex behavior of the model. Domany and Kinzel showed that this model presents a continuous transition between 
a "frozen" phase where all the sites are in the death state (this is an absorbing state from which the system cannot 
scape once it is reached) and an "active" phase with a constant average fraction of sites in the active state. Some years 
after Martins et al showecH, using the spreading of damage technique, that the active phase is divided two regions: 
one in which the system presents sensitivity to the initial conditions and one in which it does not. Such regions were 
^ , interpreted as dynamical phases and the transition between them was shown to be a continuous one. The sensibility 
'"^ ' to initial conditions is studied by looking at the time evolution of the Hamming distance bet'Vjffien two replicas of the 
^ system (i.e., two identical systems with different initial configurations) under correlated noiseu. Recently, Hinrichsen 
O et al showed that the active phase is actually subdivided into three dynamical phases which depend on the correlations 
between the different types of stochastic noise used for each replica. The phases are the following: (z) one in which 
the damage spreads for all possible types of correlations, which coincides with the so called chaotic phase founded by 
Martins et al; (ii) one in which the damage spreads for some type of noise but it does not others and [iii) one where 
, damage heals always. _ _ 

d The.jjiiiase diagram of this model was calculated using different techniques, such as mean fieldQ □ and Monte 



Carle 

In this work we introduce a generalization of the above model, in which the state of any site depends on the state 
of all the rest of the sites at the previous time through a distance dependent transition probability. In this case the 
interaction between a pair of sites separated by a distance r decays with a power law of the type 1/r" (a > 0). In 
other words, the interactions are long ranged in space. Microscopic long-range interactions appear in different physical 
contexts (see Ref.ll3 and references therein) but, until now, most of the theoretical analysis has been concentrated 
into the equilibrium properties of such kind of systems. Hence, it is of interest to analize systematically the influence 
of long range interactions in a dynamical model, such as a cellular automaton. 
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II. THE MODEL: LONG-RANGE CELLULAR AUTOMATON 



Let us consider a stochastic one-dimensional lattice with 1 sites. The configurations of the system are described 
in terms of a set of occupation variables Si = 0,1 with i = 0, 1, . . . , A^. Let S* the state of the site i at the discrete time 
t. We denote by S* = {S^} the state of the complete system at time t. We consider periodic boundary conditions, so 
that 

5*+jv+i = ^* * = 0,1,...,7V (1) 

The system evolves in time according to a spatially global rule, but without memory effects, that means, the state 
of the site i at time t depends on the state of all the rest of the sites 5*^^ at the previous time t — 1. However, the 
influence of a site j at time t — 1 over the site i at time t decays with the distance between sites with a power law 
\i — jl" with a > 0. Before establishing explicitly the evolution rule we introduce an auxiliary function p'(S'3|5'i, S2) 
of three binary variables , 6*2 and S3 , such that 

p'{o\Si,S2)+p'{i\Si,S2) = i ys,,S2 (2) 

and 

p'(l|0,0) = 

p'(l|l,0)=p'(l|0,l)=pi (3) 
p'(l|l,l)=P2 

with 



0<Pi,P2<l- (4) 

Notice that p'{S3\Si, S2) has the properties of a conditional probability for the variable S3 given the values of Si and 
S'2. Now, we consider a parallel dynamics, so that all the sites are updated simultaneously and independent of the 
rest of the sites at the same time. Then, the transition probability W{S*\S*~^) from state S*~^ to 5* is given by 

N 

WiS'\S'-')^l[w{Sl\S'~') (5) 

1=0 

where w(5'*|S'*^^) is the probability that the site i, at time t, has the value Sj given that, at time t — 1, the system 
was at S*~^. 

The model is then defined by the expression: 

W{SI\S'-') = "° P^J^d^^^n^^^+n) (g) 

where N is assumed an even number for simplicity and a > 0. Note that the interactions are symmetric on both 
sides of site i and that, together with the periodic boundary conditions (|l|) every interaction is counted just once on a 
ring of length + 1, as schematized in Fig.l. The sum J2n=^i 7^ ii^s^t appears in thedenominator of Eq.(|^) ensures 
the correct normalization of the conditional probability ti;(S'* jS**"^). It is also relatodtil to some function N*{a) that 
provides the correct scaling of the thermodynamic functions of magnetic modelsEHI for a < 1, where these kind of 
interactions are non integrable. 

From Eqs.(||)-(^) is easy to verify the following properties: 

< w{Sl\S'-^) < I (7) 

J2 wis,\s'-^) = i ys'-^ (8) 

Si=0,l 

and from Eq.(||): 

Y,w{s^s*-^)^i ys*-^ (9) 
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In the limit a — > cxd we have that 

hm wiSflS'-') = j/iSllSlzlSl-l) (10) 

a — >oc 

and the model reduces to the DK cellular automatorfl. For finite values of a the transition probability for the site 
i depends on all the rest of the sites. For small values of a the interaction strength tends to become independent 
of a. This is exactly the case for a = where the system lost its dimensionality. In the limit N —^ oo this model 
(a — 0) can also be thought as an infinite dimensional one. In magnetic systems the mean field approximation 
becomes exact in the limit of infinite dimensionality. Moreover, in a recent study of the one-dimensional Ising model 
with long-range interactions of the type l/r" it has been shown that the mean field behaviour is exact, not only in tha 
infinite dimensional case a = (which corresponds to the Curie- Weiss model) but in the full range of interactionsEiJ 
< a < 1. This result suggests that mean field behaviour could be an universal property of systems with long-range 
interactions, from which the infinite dimensional case a = is just a particular case. Therefore, it is interesting to 
investigate whether or not this property is shared by dynamical models with long-range interactions, like the present 
one. 

In a very broad sense (and perhaps a rather vague one) a mean field approximation is one in which fiuctiuaiions are 
neglected. Different systematic mean field approaches has been proposed for the DK cellular automatonB'ElLl. They 
all coincide (at least qualitatively) in the prediction for the shape of the frozen-active transition line, but they differ 
in the prediction of the damage spreading transition line. Hence, it is of interest to analize the behaviour of the 
long-range cellular automaton for small values of a, where some kind of mean-field type behaviour can be expected, 
by analogy with the magnetic models. 

Let Pt{S) be the probability of state S at time t. The time evolution of this function is given by the equation: 

p,+i(5) = ^iy(5|5')^t(5')- (11) 

S' 

We are mainly interested in the asymptotic behaviour of this model for t — > oo. Since this behaviour is not expected 
to depend on the initial state (except for some very special cases, for instance, Sj = Vi, where the solution is trivial), 
we consider, without loss of generality, the following initial state: 



^ 1 

Po{S)^l[-{Ss„o + Ss,,i) (12) 



i=0 



where 6si,Sj is a Kronecker delta function. 

We denote by {f{S))^ the mean value of a state function f{S) at time t: 

(/(5)),EE^/(5)P,(5) (13) 
s 

The mean activity is then defined as 

1 / ^ \ 1 ^ 



\i=0 

where Pt{Si) is the marginal probability for Sf. 

PtiS,)^ J2 PtiS). (15) 

Since we are considering an homogeneous system Pt{Si) is independent of the site i. Then, defining Pf(l) = Pt{Si = 
1) we have that 

ajv(t)-P(l) (16) 

This kind of system is expected to present two type of homogeneous stationary states for t —> oo: one with 5*^ = 
Vi, which defines the so-called "frozen" or "dry" phase, and one with a finite fraction of active sites Si = 1, which 
defines the "active" phase. The order parameter for this phase transition is given by 
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qn = lim aN(t). (17) 

t — *oo 

From Eqs.(||), (^ and ( p^ ) the time evolution of aNit) is given by the equation: 

aN{t + l)^{w{S,^l\S% (18) 

for an arbitrary site i. Let Pt(Si, Sj) the joint marginal probability for Si and Sj. From Eqs.(||), (^ and (|l^) we have 
that 



-9 ^^^1 

aN{t + 1) = 2pi a^(t) + ^^/,^J_ E 1) (19) 



where 

P,^{1,1) = Pt{S,^n = l,5.+„ = 1) 

and we have used that 

Pt(l,0) = P,(0,1) = Pt(l)-Pi(l,l) 

for all pair of sites i ^ j- 



III. THE FROZEN-ACTIVE PHASE TRANSITION 



Before proceeding with the mean-field and Monte Carlo calculations we can obtain a few analytical properties of 
the phase diagram in the (pi,P2) space, regarding the frozen-active phase transition. 

First of all, notice that, from Eqs.(||) and (^, the absorbing state ajv = {Si = Vi) is always a fixed point of the 
dynamics for all values of a and N . That means, a-^(t) — implies that af^{t -I- 1) = 0. This state is expected to be 
an attractor of the dynamics for low values of pi, as occurs in the DK automaton [a — oo). 

On the other hand, for the particular case P2 = 1 the fully occupied state a^^ — \ (Si = \ Vi) is also a fixed point 
of the dynamics for all values of a and N . If we now perform a "particle-hole" transformation, by defining a new 
set of variables 5* ■ = 1 — Si 'ii, we obtain, from Eqs.(||)-(^ and Eq.(^), an equivalent system with parameters P2 = 1 
and p'l = 1 — pi. We also see that = 1 — oat. Therefore, if the system presents a frozen-active phase transition 
for P2 = 1 then the critical value of pi has to be pi = p'l = 1/2. Moreover, in the active phase ajv = 1 and the 
phase transition will be a discontinuous one. These results will be verified by both the mean-field and Monte Carlo 
calculations. 



A. Mean-field solution 



From Eq.(|l9|) we see that the time evolution of ajv(t) depends on higher order correlation functions between sites. 
A simple dynamical mean-field approximation can be obtained by neglecting such correlationa3, that is, 

Pt{S,, S,) = Pt{S,)Pt{Sj) V I, J (20) 

which implies that 

Pl\l^) = mi)f = al{t) (21) 

Replacing into Eq.(p^ we find that 

aN{t + l) = 2piaN{t) + (p2-2pi)a%{t) (22) 

This equation is independent of a and N, and it is the same already encountered (with the same approximation) 
for the DK automatorfl. The fixed point equation ajv(t + 1) = aN{t) = apf derived from Eq.(p2|) gives the critical line 
p1 = 1/2 \/p2 with 

UN = (23) 

2pi-p2 
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in the active phase pi > pi- 

We can now develop a systematic improvement of this mean-field approximation, by neglecting the pair correlations 
only between sites located at a distance greater than K of site i, for some integer value K, that is, 

Pf(l,l) = [Pt(l)]' ioTn>K (24) 



with 1 <K < N/2. Then, the Eq.(|22|) can be rewritten as 



a^{t + 1) = 2p, a^it) + {p, - 2p,)al{t) il - J + E ^^ti^. 1) (25) 



For large lattice sizes ^ 1 the sum diverges for < a < 1, while 

N/2 



C(a) for a > 1, (26) 



n—l 

where C(a) is the Riemann Zeta function. Hence, for any finite value of if, Eq.(p5|) reduces to Eq.(p2|) in the limit 
— > oo, for < a < 1. This fact suggests that the mean field solution becomes exact for an infinite lattice when 
< a < 1, while a departure from such solution should be expected for a > 1, a behaviour that has already been 
observed in the ferromagnetic Ising model with long-range interactions, as was mentioned in the preceeding section. 
In fact, the Monte Carlo simulations of the next subsection confirm such conjecture. 



B. Monte Carlo results 



We performed a Monte Carlo simulation for several system sizes, N = 100, 200, 400 and 800, and for different 
values of a. It is worth mentioning that the computational complexity of the algorithm implied by the long-ranged 
transition probability (^) is 0{N^). 

We calculated the transition line between the frozen and active phases, founding a continuous transition for all 
values of a in the full plane {pi,P2), except for p2 = 1, where the numerical results verify the prediction of the 
beginning of this section. That is, for p2 = 1, the transition point always occurs at pi = 1/2 for all values of N and 
a, and the transition is a discontinuous one at that point. The typical shapes of the transition lines between the 
frozen (left part of the diagram) and active (right part of the diagram) phases are shown in Fig. 2 for N = 800 and 
different values of a. The curves show all the same qualitative shape as the corresponding one for the short-range 
model a = oo (DK model). We can observe that, for a running from a = cxd to a = 1 the curves move smoothly to 
the left of the diagram. All the curves with a between and 1 are indistinguishable inside the error intervals of the 
calculation (the error bars in the pi direction are approximately of the order of the symbol size; there are no error 
bars in the P2 direction). This behaviour was verified for all the values of A^ used in the simulation. 

For each checked value of a > 1 the corresponding curves show a fast convergence into an asymptotic one for 
increasingly higher values of A^. On the other hand, for < a < 1 the curves show a slow convergence towards the 
Pi = 1/2 vertical line, as depicted in Fig. 3 for a = 0. In fact, we verified for several values of P2 the asymptotic 
behaviour — 1/2 ~ A{p2) N~^^^ for high values of A^, where ^(^2) is some continuous function oip2 with A{1) = 0. 
The Fig. 4 illustrates the typical behaviour of pi as a function of a, for several values of A^ and a particular value of 
P2. The circles in that figure correspond to a numerical extrapolation for N 00. 

In Fig. 5 we show a numerical extrapolation for A^ — > 00 of the full critical lines for different values of a. 

Finally, in Fig. 6 we show a numerical extrapolation for A^ ^ 00 of the mean activity a — Zimjv-tooaiv pi curves 
(symbols) for three different values of p2 and < a < 1 (the activity curves for different values of a in this range are 
indistinguishable). The solid lines are the mean field prediction of Eg. ([2^). The excellent agreement between these 
results confirm the conjecture of the preceeding subsection. 



IV. SPREADING OF DAMAGE 



To study the spreading of damage in this model we have to create a replica of the system with an initial "damage" , 
by flipping randomly a fraction p of sitescl. Then we let both replicas to evolve under correlated noise and analize the 
long time evolution of the normalized Hamming distance between replicas. If the Hamming distance evolves to zero 
it is said that the damage "heals" ; otherwise, the damage "spreads" . 
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In a recent work, HinrichsenO et al have shown that the damage can spread in the DK automaton for different 
kinds of correlated noise. In this work we let the replicas ev(jJ-Me,uiider the identical realizations of the stochastic 
noise, which was the most studied case in the DK automatorO'BQ&El. Thia, corresponds, according to Hinrichsen et al 
definition, to maximally correlated noise, at least for the DK automatond {a = oo). It is expected this property to 
be valid for all a, as long as in this case a single random number is generated for every site in both replicas. Other 
types of correlated noise are very difficult to simulate in the long-range model, because they involve, in principle, the 
generation of 0{2^) random numbers for every site at every time step. 

We performed this Monte Carlo simulation for iV = 100, 200 and 400 founding that, for < a < 1, the damage 
spreading is suppressed for all values of {pi,p2)- For a > 1 a small region where the damage spreads appears in the 
lower right corner of the diagram. The transition line moves continuously towards the DK one for a = oo. 

In Fig. 7 we show the transition lines between the healed and the "chaotic" regions, for different values of a > 1 
and N = 400. 



V. CONCLUSIONS 



We introduced a one-dimensional stochastic cellular automaton which allows to analyze different ranges of interac- 
tions through a single parameter a. We verified that this model reproduces, for high values of a most of the known 
results of the DK model (a = oo). 

We found two different regimes, according to the range of values of a. 

For a > 1 the dynamical phase diagram shows the same qualitative aspects of the one corresponding to the DK 
model, that is, three different regions: a frozen and an active phase, with the last one divided in two regions, one 
where the damage spreads ("chaotic" phase) and one where the damage heals, when the replicas evolve under the 
same noise. The transition lines between these regions continuously vary when a decreases from a = oo to a = 1. 

For < a < 1 the phase diagram becomes independent of a, even for finite N. No damage spreading were found in 
this case, at least for replicas evolving under the same noise. Regarding the frozen-active phase transition, we presented 
a strong evidence that the mean field theory is exact in theJimit N oo. The fact thaLthe same behaviour appears 
in the thermodynamics of the one-dimensional Ising modelllii with similar interactional^ may suggests the existence 
of some universal relationship between mean field theories and very long ranged interactions (0 < a < d, where d is 
the space dimension). 

It is worth noting that mean field caJeulations of the Hamming distance dynamics between replicas evolving under 
two different types of correlated noiscB-El predicts a spreading of damage transition in the DK model. It is a general 
heuristic argument (based on the equilibrium behaviour of magnetic systems) that mean field approximations of short 
range models (a — oo) show the same behaviour as the corresponding fully connected models (a — 0). We have 
shown that this is not true for the spreading of damage transition in the DK model. This variance between the two 
phenomena in the cellular automata is telling us that the spreading of damage transition is of different nature, at 
least in some of its properties, than other cooperative phenomena like the frozen-active or the ferromagnetic phase 
ones. l-. 

Finally, there is recent conjecture of Grassbergeia which states that all continuous transitions from an "absorbing" 
state to an "active" one with a single order parameter are in the universality class of directed percolation, provided 
that some conditions are fulfilled, among others, short-range interactions in space and time. Both transitions in 
the present model fall into this category, except for the range of the interactions for small values of a. Hence, having 
a model which allows to interpolate between long and short range regimes, and for which the long range behaviour is 
known exactly may be a good candidate for testing the above conjecture. 

I am indebted to Constantino Tsallis for many fruitful suggestions and discussions about this problem. Useful dis- 
cussions with M. L. Martins, T. J. Penna and M. de Oliveira are acknowledged. I also acknowledge warm hospitality 
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CONICET (Argentina), Consejo Provincial de Investigaciones Cientfficas y Tecnologicas CONICOR (Cordoba, Ar- 
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There is not a simple mapping for finite a between this model and the equilibrium thermodynamics of a magnetic model, as 



occurs m 



the DK modelM. 



FIG. 1. Schematic diagram of interactions for the long-range cellular automaton in a ring of A'^ + 1 sites with periodic 
boundary conditions. The probability of site i to be in a given state at time t depends sym metrically on all the rest of the 
sites along the ring (at time t — 1), down to the site i — N/2 and up to the site i + N/2 {N is assumed an even number). 



FIG. 2. Monte Carlo calculation of the transition lines between the frozen (left part of the diagram) and the active (right 
part) phases for A'^ — 800 and different values of a. The lines are a guide to the eye. 



FIG. 3. Monte Carlo calculation of the transition lines between the frozen (left part of the diagram) and the active (right 
part) phases for a — and different values of N. The lines are a guide to the eye. 



FIG. 4. Monte Carlo calculation of the critical value pi as a function of a for p2 — 0.5 and different values of N. The circles 
correspond to a numerical extrapolation for N — » oo. The error bars, when not shown, are smaller or equal than the symbol 
size. The lines are a guide to the eye. 



FIG. 5. Numerical extrapolation for A'^ ^ oo of the critical lines for different values of a. 



FIG. 6. Numerical extrapolation for N <x of the mean activity a (symbols) compared with the mean field prediction 
(solid lines) for < q < 1 and different values of p2 . 



FIG. 7. Transition lines between the healed (left part of the diagram) and the "chaotic" (right part) phases, for different 
values of a > 1 and A'' = 400. 
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